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حکیده 


یکی از مهم‌ترین چالش‌های محیطی در سالیان اخیر در استان خراسان رضوی پدیده گرد و غبار است. بر 
این اساس پیش‌بینی فضایی- زمانی تعداد روزهای گرد و غباری سالیانه با استفاده از روش کریجینگ و با 
کمک نرم‌افزار 15 مورد بررسی قرار گرفت که در آن سرعت باد ۱۵ متر در ثانیه و بیشتر و دید افقی زیر 
۰ متر از نظر سازمان هواشناسی جهانی به عنوان روز گرد و غباری در نظر گرفته شد. پس از حذف 
همپوشانی داده‌ها از دو ماتریس فضایی- زمانی دید افقی و سرعت باد. آرایه 12212 ٩۳‏ به صورت ترکیبی 
از ماتریس و بردار در کلاس ۳ 91۳ ساخته شد. سپس تمامی مدل‌های تفکیک‌پذیر و غیر 
تفکیک‌پذی به مدل تجربی داده‌ها برازش داده شدند که تغییرنگار متریک با کمترین میانگین مربعات خطا 
به عنوان بهترین مدل برای پیش‌بینی تعداد روزهای گرد و غباری سالیانه انتخاب شد. خروجی مدل نشان 
داد که داده‌ها تا ۵ سال دارای وابستگی فضایی- زمانی هستند, لذا می‌توان تا سال ۲۰۲۲ تعداد روزهای 
گرد و غباری را برآورد نمود. حدود اطمینان تعداد روزهای گرد و غباری در سطح ۸4۵ نشان داد در سال 
۵۹ ایستگاه سبزوار با ۵1 روز بیشترین و ایستگاه گناباد با ۲۳ روز کمترین تعداد روز گرد و غباری را 
خواهند داشت. همچنین پیش‌بینی‌ها نشان می‌دهد تعداد روزهای گرد و غباری در استان خراسان رضوی 
رود کا و وق ال ۰۲ ی وا که تلا اف نش شارت 

کلیدواژه‌ها: تغیبرنگار فضایی- زمانی» پیش‌بینی. روزهای گرد و غباری. استان خراسان رضوی. تحلیل‌های 
فضایی - زمانی. 
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۳21 جغرافیا و مخاطرات محیطی شماره سی و دوم 


۱- مقدمه 

ازجمله مهم‌ترین پدیده‌های هواشناسی در نواحی استان خراسان رضوی می‌توان به طوفان شن و خاک" با سرعت 
باد ۱۵ متر بر ثانیه و بیشتر و طوفان گرد و غبار یا ریزگرد" با دید افقی کمتر از ۱۰۰۰ متر اشاره کرد (سازمان 
هواشناسی جهانی؛ ۲۰۱۵). به این دو پدیده به اصطلاح طوفان گرد و غبار " گفته می‌شود (سازمان هواشناسی جهانی؛ 
۱ سازمان جهانی بهداشت برآورد نموده است که سالانه ۵۰۰ هزار نفر بر اثر مواجهه با ذرات معلق هواپرد 
موجود در هوای آزاد. دچار مرگ زودرس می‌شوند. بر اساس این برآورده هزینه سالیانه صرف شده برای بخش 
سلامتی و بهداشت ناشی از آلودگی هوا در اتریش. فرانسه و سوئیس حدود ۳۰ میلیارد پوند بوده است که اين میزان 
معادل 1 درصد از کل مرگ و میرهای این کشورها است (1.2003 5ع[0)). 

این وضعیت در سال‌های اخیر در استان خراسان رضوی شرایط مشابهی داشته است. تحلیل فضایی - زمانی 
کدهای طوفان‌های گرد و غبار نشان می‌دهد که میانگین روزهای گرد و غباری در ایستگاه‌های منتخب استان خراسان 
رضوی طی دوره آماری ۱۹۸۷ تا ۲۰۰7 برابر با ۲۱ روز بوده است که مشهد با ۳۸ روز بیشترین» سبزوار و 
تربت‌حیدریه با ۱۶ و ۱۰ روز کمترین تعداد روزهای گرد و غباری را داشته‌اند (فرج زاده و علیزاده۱۳۹۰). همچنین 
تحلیل آماری- همدیدی طوفان‌های گرد و غبار در فاصله زمانی ۱۹۹۳ تا ۲۰۰۵. نشان می‌دهد تعداد طوفان‌های گرد 
و غبار از شمال به جنوب استان حراسان رضوی افزوده می‌شود. (لشکری و دیگران۱۳۸۷) با این وجود بررسی‌ها در 
ایستگاه سینوپتیک سبزوار در طی دوره آماری ۷ تا ۱۳۸۷ نشان می‌دهد بیشترین احتمال وقوع گرد و غبار در 
ماه‌های اردیبهشت و خرداد اتفاق افتاده است و بیش از ۷۱ درصد در بعدازظهرها روی می‌دهد. همچنین وقوع 
روزهای همراه با گرد و غبار در شهرستان سبزوان روند افزایشی داشته است (مهرشاهی و نکونام ۱۳۸۸)» و تعداد 
روزهای گرد و غباری آن نیز به‌طور متوسط در سال به ۱۵ تا ۱۷ روز می‌رسد که فصل بهار شایع‌ترین زمان بروز این 
پدیده بوده است (امیدوار و نکونام» ۱۳۹۰). همچنین مطالعات در دوره آماری ۲۰۰۶ تا ۲۰۱۳ از نظر پراکندگی فضایی 
در استان خراسان رضوی نیز نشان می‌دهد. بیش‌ترین روزهای گرد و غباری مربوط به نواحی شمال شرق سرخس به 
دلیل نزدیکی با بیابان ترکمنستان» شرق و غرب و جنوب استان خراسان رضوی بوده است که کاهش پوشش گیاهی 
بر تعداد روزهای گرد و غباری تأثیر داشته است؛ به‌طوری که در سال ۲۰۰۸ بیش‌ترین روزهای گرد و غباری دیده 
شده است که این سال خشک‌ترین سال دوره آماری بوده است (پورهاشمی سمیرا و همکاران ۱۳۹۶). 

مطالعات پایش ماهواره‌ای و تحلیل همدید پدیده گرد و غبار در کلان شهر مشهد. طی دوره آماری ۲۰۰۹ تا ۲۰۱۳ 
برای تعیین مناطق منشاً گرد و غبار با کمک مدل روش پسگرد در ۲۶ ساعت قبل از وقوع در روزهای آماری, بسته 
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سال هشتم پیش بینی تعداد روزهای گرد و غباری سالیانه در ... ۱۵ 


به مسیر گرد و غباری نشان داد در طی دوره گرم سال آسیای مرکزی. بیابان‌های شرق ترکمنستان, افغانستان و اراضی 
تغییریافته دشت‌های شمال شرق کشور و در طی دوره سرد سال طوفان‌های گرد و غباری جهت خاصی ندارند. 
۱ 

استفاده از مدل شبکه عصبی مصنوعی در بررسی روابط بین عناصر اقلیمی و نوسانات آلودگی هوا در کلان شهر 
مشهد (کاخکی مهنه» ۱۳۹۲ پیش‌بینی ارزش گذاری و ترسیم نقشه آلودگی هوای شهر مشهد (ضیائی و همکاران؛ 
۰ ) بررسی آماری سینوپتیکی آلودگی هوای مشهد (جهانشیری» ۱۳۸۹ ارزش گذاری ویژگی‌های مختلف 
آلودگی هوا در مشهد (قربانی و همکاران» ۱۳۸۸) و بررسی سینوپتیک اینورژن در مشهد با استفاده از تحلیل عاملی 
(علیجانی و نجفی نیک ۱۳۸۸) مواردی از تحقیقات انجام شده در این منطقه می‌باشد که هیچ کدام از اين تحقیقات بر 
مدل‌های پیش‌بینی تاکید نداشه‌ند؛ لذا با توجه به اهمیت موضوع و فراوانی وقوع پدیده گرد و غبار در استان خراسان 
رضوی هدف اصلی از این پژوهش پیش‌بینی تعداد روزهای گرد و غباری سالیانه با کمک مدل فضایی- زمانی است 
تا بتوان نقاط بحرانی را در زمان‌ها و مکان‌های مورد دلخواه پیش‌بینی و از خسارت وارده به این مناطق جلوگیری کرد 
که با فرض معلوم بودن پارامترهای مدل» شامل میانگین؛ واریانس و در نظر گرفتن تغییرات توأم فضایی- زمانی 


داده‌ها؛ پیش‌بینی تعداد روزهای گرد و غباری سالیانه در سال‌های آینده مورد بررسی قرار می‌گيرد. 


۲- مواد و روش‌ها 

۱-۲- منطقه مورد مطالعه 

محدوده مورد مطالعه به ترتیب بین ۵0 تا ٩۱‏ و ۳۶ تا ۳۸ درجه طول و عرض جغرافیایی در شمال شرق کشور 
واقع شده است (شکل ) موقعیت مراکز شهرستان‌ها و ایستگاه‌های مورد مطالعه در استان خراسان رضوی و 
ایستگاه‌های مجاور در شکل دو آمده است. 
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شکل ۱- منطقه مورد مطالعه شکل ۲- موقعیت ایستگاه‌های مورد مطالعه 
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شکل شماره ۳ میانگین روزهای گرد و غباری در طول دوره آماری نشان می‌دهد که ایستگاه فریمان با ۲۵ روز و 
خواف با ۲۳ روز. گلمکان با ۲۱ روز و مشهد با ۱٩‏ روز دارای بیشترین و ایستگاه‌های سبزوان نیشابور تربت‌حیدریه, 
بردسکن. گناباد تربت جام و کاشمر کمترین تعداد روزهای ق و غباری سالبانه را در طول ففدم آمادع داشته‌اند. 


میانگین تعداد روزهای گرد وغباری درایستگاههای سینوپتیک استان خراسان رضوی 
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شکل ۳- میانگین تعداد روزهای گرد و غباری سالیانه در طی دوره آماری ۱۹۸۷/۱۱/۰۱ لغایت ۲۰۱۳/۱۲/۳۱ 


شکل ۶- تعداد وزهای گرد و غباری سالیانه در استان خراسان رضوی را نشان می‌دهد که علیرغم نوسانات 
شدید. روند کلی آن در سال‌های اخیر در حال افزایش است. نتایج برازش مدل رگرسیون خطی نشان داد که این 


روند معنی دار (آماره ۳ برابر ۶6/۱۲ و مقدار برابر ۳/۳۱۱۶>۱۰۳۲) بوده است. 


میانگین تعداد روزهای گرد و غباری سالیانه در استان خراسان رضوی 
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شکل - روند روزهای گرد و غباری در استان خراسان رضوی 


سال هشتم پیش بینی تعداد روزهای گرد و غباری سالیانه در ... ۱۷ 


۲-۲- روش پژوهش 

در این تحقیق از آمار سرعت باد ۱۵ متر و بیشتر و دید افقی کمتر از ۱۰۰۰ متر به ترتیب به عنوان 5۸ و ل1«0 
دردوره آماری ۱۹۸۷/۱/۱ تا ۲۰۱۷/۳/۳۱ به صورت روزانه استفاده شده است که با در نظر گرفتن ۱۶ ایستگاه در 
منطقه مورد مطالعه پس از استاندارد سازی داده‌هاء در میدان تصادفی ()2 به صورت زیر تجزیه شده است 
(1993 .۱۲ ,065516)) 


۲ بر( ر۵)5 + ارو > (ارو)7 


جهت مدل سازی, پردازش و برآورد داده‌ها از دو فرمت ٩1۳‏ "و 5۲۳۲۳ ۲ (2013 .ظ۴ ,۵۵92ع۳) در طی 
سال‌های پیش رو استفاده شد که هر کدام از این دو فرمت ترکیبی از ماتریس و بردار است بدین سان داده‌های فضایی 
و زمانی در قالب جداول جداگانه به صورت 2,۷,۲ ساخته شدند سپس در فرمت‌های طولانی با مشخص کردن هر 
رکورد» ترکیب زمانی و فضایی داده‌ها" به دست آمد (۵۵9۳02,.2012ع۳) 

با توجه به منطقه مطالعه و هدف پژوهش با حذف مقادیر داده‌های مفقود شده* شکل داده‌ها تحت عنوان آرایه 


مق ٩0‏ با توجه به رابطه زیر ساخته شد: 


رمشاک زو رف لیاوا 
در آن اندیس ژبه شماره مکان (موقعیت مشاهده) و اندیس ژ به شماره زمان اشاره دارد. 


سپس بردار 26,۷ داده‌ها به صورت سیستم تصویر مرکاتور معکوس جهانی " تهیه شد. فاصله قطری" مستطیل 
محاطی محدوده مورد مطالعه از رابطه: 


[2/3۴ 


محاسبه شد در نهایت برای تغییر نگار تجربی زمانی " تا ۱۰ سال و برای تغییر نگار تجربی فضایی" تا 1۵۰ 


کیلومتر در کلاس 8۳ محاسبه گردید. سپس با توجه به تابع کواریانس فضایی-زمانی ": 
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۱۸ جغرافیا و مخاطرات محیطی شماره سی و دوم 
(ا)ع(وت)ی) < (۷۵,وال)) 


که در آن: 

(ود)ی) تغییر نگار صرفاً مکانی 

و(م)) تغییر نگار صرفاً زمانی است. 

تابع هم تغییرنگار هر کدام از تغییرنگارها به تفکیک به دست آمد و داده‌ها به دو تغیبرنگار فضایی و زمانی به 
صورت مجزا برازش داده شدند (61,00.2001 069876 126) پس از آن تغییرنگار تجربی فضایی و زمانی با استفاده از 
روش کریجینگ و از رابطه زير محاسبه شد (محمدزاده :)۱۳٩۱‏ 


که در آن: 


2 ۳ 
۱ بوز ۳ 6 ۱6 , ( ر5) 7 - ( + ۲+ 5( ,»10 ت 7/9 


(مول ۷ 

که (۸۸,و0) ۷] نشان دهنده مجموعه همه زوج مشاهداتی است که فاصله آن‌ها در همسایگی بردار ء و فاصله 
زمانی آن نزدیک ,8 باشد. 

پس از ترسیم تغیبرنگار فضایی - زمانی تجربی داده‌هاء جهت برازش داده‌ها به مدل تغییرنگارضایی- زمانی 
استقلال و هتکن داده‌ها بررسی شد. برای بات کی فضایی از آزمون موران (.2010 2 :08/ونا) 6 0 بطفاع02) 
وتزای همیستکی زهاتی از را ند تست (2015 .۸۵ ,1/۵185 ۴.6 ,0ت6)) استفاده شد. جهت تعیین نوع کریجینگ, 
مانایی و همسانگردی داده‌ها نیز مورد بررسی قرار گرفت و نقشه سه بعدی تغییرنگار نیز ترسیم شد. سپس مدل‌های 
نظری گوسی. نمایی» کروی» سهمیگون. خطیء بسل و مترنبه مدل‌های تغییر تغییر نگار تجربی فضایی و زمانی 
جداگانه برازش داده شدند و در نتیجه مقادیر آستانه جزئی اثر قطعه‌ای. دامنه و مجذور مربعات خطا محاسبه شد. 
پس از آن تمامی مدل‌های نظری فضایی - زمانی تفکیک‌پذیر و تفکیک ناپذین مجموعاً ۱۶۶ مدل» شامل: ضربی» 
جمعی - ضربی متریک. متریک جمعی و متریک جمعی ساده با تغییر نگارتجربی فضایی- زمانی برازش داده شدند 
(2016,.ظ ,بعلقتن ک رمصوع0ه۳) که مدل ضربی - جمعی با داده‌هاء همگرا ند و در نتیجه اين مدل حذف 
گردید. در نهایت بهترین مدل با کمترین مجذور مربعات خطا جهت برآورد و پیش‌بینی تعداد روزهای گرد و غباری 
سالیانه انتخاب شد. در نتیجه برآورد تعداد روزهای گرد و غباری با استفاده از روابط زیر به دست آمد (محمدزاده 


۵9 


1۳( ,2)6] - (وا ,وع) ۳02 
(ره)|رمء ,۳)2)5۵ كِ "[(0) ,2)50] 
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که در آن» عبارت رر)| رمع رن5))نل مقدار مورد انتظار رم «7,)5 با شرط در اختبار داشتن 1 مقدار آن برایر با 
( رپرک) (و- ۳ رپر7,)5 قرط رو7)5 است. این بدان معناست که برآورد کریجینگ به صورت تعداد روز و 
غباری سالیانه خواهد بود. در این راستا کریجینگ برای موقعیت فاقد مشاهده (20 ,5) به صورت زیر به دست آمد 


(محمدزاده ۱ ۳۹ 


۷ ۲ 
(6(۳)۵,۵ :۱6 2 - [7)5,,00] 
12-1 1<ز 
که در آن» (:2)5::6 ضرایب کریجینگ هستنده که.در نهایت بر اساس کلاس 81۳10۲ محاسبه و تغییرنگار سه 

بعدیو رویه آن نیز ترسیم شد (2015 .1,61 ,ا۳6080) پیش‌بینی برای 7 سال آینده به صورت عددی و بر اساس 
سیستم مختصات 71 صورت پذیرفت و مقادیر پیش‌بینی شده بر روی منطقه مورد مطالعه پلات شدند. یکی از 
مهم‌ترین شاخحص‌ها جهت حصول اطمینان از برآورد احتمال روز گرد و غباری. محاسبه مقادیر خروجی در سطح 
اطمینان /٩۵‏ است که با استفاده از رابطه ٩۳(‏ .06551607): 


[(1:960:6۵0 + (ع26) ,»1:96 - (و2)۵) < ۸ 


برای تک تک ایستگاه‌ها در استان خراسان رضوی انجام شد. این موضوع حداکثر قابلیت وقوع تعداد روزهای 

گرد و غباری (کران بالای فاصله اطمینان ۹۵/) را نشان می‌دهد که می‌تواند میزان حصول اطمینان از مدل عددی را 
بدین منظور از بسته‌های نرم افزاری ]ماوع (2017 بلق 61 مقعطعع‌هعظ), عمصتام)۹۵۵( ۱۲ 21,۲۰ ۵1 رقعصعع‌هعظ) و ظ 
(2017 ,له 6 ,قصععطاع) عتفه۳ (1,2016 تعحامعل. مود (2017 .له ی بظ ,صهتنظ) و ومد عاعممن0 ۲ 
(10270.171,2013ت16 و 2016 ۷۲۰ ,10601067) برنامه نویسی لازم انجام و از مح محیط نرم افزاری 1 ( 1۵۷۵100۴60۲ 1 


1 :۲627 00۳6)) _ استفاده شده است 


۳- نتایج و بحث 

۱-۳- آزمون و برازش مدل تغییر نگار 

با برازش داده‌ها به مدل تغییرنگار فضایی- زمانی» استقلال و همبستگی داده‌ها نشان داد که داده‌های ٩15‏ در 
سطح اطمینان ۰/۰۵/ معنی دار و وابسته هستند. لذا می‌توان روزهای گرد و غباری را برای سال‌های آینده پیش‌بینی 
کرد. مشاهدات در محورهای مختصات و ۷ وا برای طوفان شن و طوفان گرد و بیانگر این است که داده‌ها دارای 


۱۳۹۰ جغرافیا و مخاطرات محیطی شماره سی و دوم 


روند خاصی در هیچ یک از مقادیر (نا,ری) 2 در جهت‌های و ۲و انیستند (شکل ۵ مقادیر ۷ با تغییررنگ 
داده‌ها نشان می‌دهد که داده‌ها دارای مانایی ذاتی بوده و می‌توان بر روی آن تحلیل فضایی- زمانی انجام داد ( طول 


جغرافیایی»۳06] زمان» ۷ تغییرات رنگ 72 مقدار عددی متغیر). 


هد که5 امنصصد که ۲حص6ااه5 


۷۵۲۱۵۵ ادتادمه آه م۵ و کی 109۳0 
1 


لفت 5 
مه 
5000 تِ 5 3 ‌ 
ٍ "۳ ۶ رب 1- 500 
۳ و تسه 
3000 * ۳ ‌ 
۳۳ كت ‌ 1 0 < 
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1000 اُ ۳ 1 
كت ۳ 3 ‌- 
ً ۳ 500 
500 ۷ 00- 
1 
۰ ۳ ۰ 2 ۰ .2 ۳۲ 
شکل ۵- بررسی مانایی در داده شکل *- بررسی همسان کُردی تغییرنگار فضایی 


شکل " رویه نیم تغییرنگار فضایی را نشان می‌دهد. با توجه به اين که نیم تغیبرنگار در تمام جهات تقریباً بر هم 
منطبق است و تفاوت معنی داری نسبت به هم ندارد؛ لذا داده‌ها همسانگردند. 

۲-۳ برآورد تغییرنگار فضایی - زمانی 

تغییر نگار تجربی فضایی - زمانی (ع)ون) تا فاصله 1۵۰ کیلومتری نشان می‌دهد که گاماء سیر صعودی داشته و 
در لگ های انتهایی» نوسانات شدیدی در آن دیده می‌شود. دامنه این تغییر نگار طولانی است و در نهایت به سقف 
نمی‌رسد. تغییر نگار تجربی زمانی (0)) با فاصله زمانی ۳۰ سال نشان می‌دهد علیرغم اين که مقادیر خروجی گام 
در لگ های میانی به هم نزدیک‌تر می‌شود در نهایت دامنه آن نیز طولانی است در هر دو تغییرنگان بزرگ‌تر بودن 


آستانه جزئی از اثر قطعه‌ای شرایط خوبی برای برازش مدل دارد. 
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شکل ۷- مدل مترن برای تغییر نگار فضایی ۸- مدل مترن برای تغییرنگار زمانی 


جهت برآورد و پیش‌بینی پارامتر نامعلوم روزهای گرد و غباری مدل‌های گوسی, کروی. خطیء بسل, مترن و 
نمایی به تغیبرنگار تجربی فضایی و تغیبرنگار تجربی زمانی و به تمامی مدل‌های تفکیک‌پذیر و غیر تفکیک‌پذیر شامل 
مدل‌های ضربی» ضربی جمعی» متریک متریک ساده و متریک جمعی؛ مجموعاً :۱ مدل برازش داده شدند (جدول 
۱( 


جدول ۱- برازش مدل‌های فضایی - زمانی تجربی با مدل‌های نظری تفکیک‌پذیر و غیر تفکیک پذیر 


مدل متریک جمعر مدل متریک ساده مدل متریک مدل ضربی مدل فضایی-زمانی 
۱,۹( ۱۹ ۱۳۳ ی( مترن+ مترن 


که مدل مترن با ضریب کاپای ؛ برای تغییر نگار حاشیه ایی فضایی (شکل ۷) و مدل مترن با قنت نس 
کاپای 0.6 برای تغییر نگار حاشیه ایی زمانی (شکل ۸) انتخاب شدند. خروجی‌ها نشان می‌دهند مدل متریک بهترین 
مدل برازش داده شده جهت برآورد سالیانه روزهای گرد و غباری است که با کمترین مجذور مربعات خطا همراه 
است. با توجه به دو مدل تجربی فضایی و تجربی زمانی. تغیبرنگار تجربی سه بعدی (شکل )٩‏ و بهترین مدل 


برازش داده شده آن نیز (شکل ۱۰) ترسیم شد. 


۱۱۲ شماره سی و دوم 


شکل -٩‏ تغییر نگار تجربی سه بعدی 


تغییرات در لگ های زمانی و فاصله افقی نشان می‌دهد. داده‌ها هم در جهت زمانی و هم در جهت فضایی؛ 
تاش ها ی کر اب کی را سس «ازفن 

۳-۳- پیش‌بینی فضایی- زمانی 

ساختار همبستگی فضایی- زمانی داده‌ها نشان می‌دهد که تغییر نگار فضایی - زمانی متریک می‌تواند طوفان گرد 
و غبار را در استان خراسان رضوی پیش‌بینی کند با این وجود در این پژوهش بازه‌ی ۵ ساله از سال ۲۰۱۸ تا سال 
۲ مد نظر قرار گرفته است که نتایجم خروجی فضایی - زمانی آن در شکل شماره ۱۱ و جدول شماره ۲ آمده 
انبست: تشیرات سای کر براکتشن فضای‌عومانی فاده‌ها فر طرل سال‌های موزد افظار اژسال ۲۰۱۸ فیله مي شوه 
بنابراین مهم‌ترین نقاطی که ارقام پیش‌بینی شده تعداد روزهای گرد و غباری سالیانه از مقادیر بالاتری برخوردارند 
عبارتند از: شمال شرق. شمال غرب و جنوب شرق استان خراسان رضوی. 

مقادیر پیش‌بینی شده سالیانه روزهای گرد و غباری در نواحی مرکزی استان خراسان رضوی در سال ۲۰۱۸ و 
۹ نشان می‌دهد که اين مناطق کمترین تعداد روزهای گرد و غباری را خواهند داشت ولی به تدریج بر میزان 
تعداد روزهای گرد و غباری آنها افزوده خواهد شد؛ به‌طوری که در سال ۲۰۲۲ اغلب نقاط استان به جز نقاط 
محدودی در شرق» غرب و جنوب شرق تعداد روزهای گرد و غباری در نواحی مرکزی و سایر نقاط افزایش 


چشمگیری خواهند داشت. (شکل 65 
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شکل ۱۱- پیش‌بینی فضایی - زمانی روزهای گرد و غباری با استفاده از مدل متریک 


۳-- برآورد واریانس خطا 

برآورد واریانس خطا برای ۵ سال مد نظر نشان می‌دهد دامنه آن از ۳۹ روز شروع شده و تا ۶۳ روز ادامه دارد که 
با توجه به دامنه محدود پیش‌بینی با محاسبه جذر واریانس خطا مقادیر بالایی را نشان نمی‌دهد. همچنین خروجی‌ها 
نشان می‌دهد که کمترین تعداد روز گرد و غباری در استان خراسان رضوی در سال ۲۰۱۵ مربوط به ایستگاه گناباد 


حداکثر با ۷ روز و بیشترین تعداد روز گرد و غباری مربوط به ایستگاه سبزوار حداکثر با ۳۵ روز است. (شکل ۱۲) 
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شکل ۱۲- حدود اطمینان تعداد روزهای گرد و غباری در سطح اطمینان ۸٩۵‏ در سال ۲۰۱۹ 


۳4 جغرافیا و مخاطرات محیطی شماره سی و دوم 


حدول۲- پیش بینی تعداد روزهای گرد و غباری و حداکثر احتمال وقوع در سطح اطمینان ۹۵ 


حد بالای | پیش‌بینی | حدبالای | پش‌بینی | حدبالای | پش‌بینی | حدبالای | پش‌یینی | حد بالای | پش‌بینی 7 شما 

اطمینان | ۲۰۲۲ | اطمینان ۲۱ اطمینان ۰ | اطمینان | ۲۰۱۹ | اطمینان | ۲۰۱۸ تفت 
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39 18 38 17 36 16 35 15 33 14 بردسکن 9 


28 7 27 7 26 7 26 7 23 6 گناباد 10 

46 24 45 23 43 ره 41 21 40 0 | تربت جام | 11 
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33 12 32 12 31 11 30 10 29 10 کاشمر | 14 

43 22 42 21 41 21 40 20 39 19 میانگین | 15 
6- نتیجه گیری 


بررسی روند میانگین تعداد روزهای گرد و غباری سالیانه در استان خراسان رضوی نشان می‌دهد که از سال 
۷ با میانگین ۱۲ روز به ۲۱ روز در سال ۲۰۱ رسیده است و از سال ۱۹۹۳ به بعد. به شدت سیر صعودی به 
خود گرفته است. به‌طور کلی پراکنش لگ های تغییرنگار تجربی نشان می‌دهد که داده‌ها هم در جهت زمانی و هم 
در جهت فضایی دامنه طولانی دارند به همین دلیل تغییرنگار فضایی- زمانی به سقف نمی‌رسد. این بدین معناست 
که تعداد روزهای گرد و غباری در استان خراسان رضوی به شدت در حال افزایش است. برازش مدل‌های مختلف 
تغییر نگار فضایی- زمانی نشان می‌دهد که مدل متریک با کمترین میانگین مجموع مربعات می‌تواند تعداد روزهای 
گرد و غبار سالیانه را پیش‌بینی کند 

تغییرات آستانه جزئی نشان داد که پیش‌بینی برای ۵ سال آینده می‌تواند انجام شود. پراکنش فضایی- زمانی تعداد 
روزهای گرد و غباری نشان می‌دهد تغییرات در بعد فضا در سال ۲۰۱۸ و ۲۰۱۹ کمتر دیده می‌شود و به نقاط 
محدودی از شمال شرق. شمال غرب و جنوب شرق استان خراسان رضوی خلاصه می‌گردد؛ اما تغییرات زمانی تا 


شسال ۲۰۲۲ نشان مي دهد اغلت تقاط اسان و تواععی مرکزی افزایش خشمگری شواهند داشت: وونل تعیی ان 
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فضایی- زمانی نشان می‌دهد در سال ۲۰۱۹ ایستگاه‌های قوچان, مشهد. سرخس, سبزوار خواف» تربت جام و 
ترا فا ۲۱۸۵۲۱۱۳۱۳۵۱۲۲۱۲ سای روز کرقی یارس هشن فنوانت که حتاف قانلرت اعسمام در 
سطح /٩۵‏ در این سال به ترتیب در همان ایستگاه‌ها به 1۸۱۵76۷۵۱6۱ روز می‌رسد که بیشترین تعداد وقوع 
روز گرد و غبار مربوط به ایستگاه سبزوار با ۵7 روز و پس از آن سرخس با 4۷ روز است. در اين سال ایستگاه‌های 
گناباد با ۷ کاشمر با ۱۰ و درگز با ۱۳ روز با کمترین تعداد روز گرد و غباری پیش‌بینی شده‌اند که حداکثر احتمال 
وفوع گرد و غبار در این ایستگاه‌ها در سطح اطمینان ۸۹۵ به ترتیب به ۰۲٩‏ ۳۰ و ۳۳ روز می‌رسد. 

روند پیش‌بینی فضایی_زمانی تعداد روزهای گرد و غباری سال ۲۰۲۲ نشان می‌دهد که ایستگاه گناباد با ۷ روز 
تعداد روز گرد و غباری چندانی را نشان نمی‌دهد و بهترین شاخص کیفیت هوا را دارد که حداکثر احتمال وقوع آن 
تا ۲۸ روز می‌رسد. همچنین حداکثر احتمال وقوع تعداد روزهای گرد و غباری در ایستگاه‌های قوچان؛ مشهد. 
سرخس, سبزوار, خواف» تربت جام و فریمان به ترتیب با تعداد,۵۱۶۵۵ ۵۱6۵۹۵۲۱ روز پیش‌بینی شده‌اند که 


می‌تواند در برنامه ریزی منطقه‌ای و ناحیه‌ای در این استان مد نظر قرار گیرد. 


کتابنامه 

آمار سازمان هواشناسی کشور. سرعت باد و دید افقی از ۱۹۸۷/۱/۱ لغایت ۲۰۱7/۱۲/۳۱ 

امیدوا کمال؛ نکونام. زری؛ ۱۳۹۰. کاربرد گلباد و گل غبار در تحلیل پدیده گرد و خاک و تعیین رژیم فصلی بادهای 
همراه با این پدیده (مطالعه موردی: شهر سبزوار). پژوهش‌های جغرافیای طبیعی. ۰۷۲ ۸۵-۱۰۶ 

پورهاشمیء سمیرا؛ بروغنی. مهدی؛ زنگنه اسدی: محمدعلی؛ امیراحمدی, ابوالقاسم؛ ۱۳۹۶. تحلیل ارتباط پوشش 
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